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Abstract—The dependence between the processes of formation of 
classes of primes in the generalized hypothesis of Artin is studied on 
the basis of the theory of randomized algorithms in algebraic and 
analytic number theories. It is proved that the construction of 
randomized algorithms for estimating Artin's constants is the basis for 
the construction of computer models of classes of primes in 
accordance with the generalized Artin's hypothesis. It has been 
established that modern methods of analytic number theory do not 
allow one to obtain estimates for the generalized Artin constants. 
Methods for calculating Artin's constants for arbitrary sets of natural 
numbers are developed and the convergence of estimates of constants 
in probability to limit values is established. The foundations of the 
number-theoretic analysis of Artin's constants and related classes are 
created. The dynamics of changes in the laws of distribution of prime 
numbers depending on the magnitude of the number of prime divisors 
of the previous number has been studied. A hypothesis is formulated 
about the limiting properties of the distribution of primitive roots of 
prime numbers. 
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I. INTRODUCTION 


The solution of many problems in various areas of applied mathematics 
depends on the solution of a significant number of problems in pure 
mathematics that have not yet been solved. Artin's primitive roots conjecture is 
one of these fundamental mathematical problems. For almost a century, its 
solution has not been brought to such a level that would allow improving 
methods for solving the discrete logarithm problem, modern cryptography 
algorithms, models of algebraic dynamical systems, and developing methods 
for constructing more advanced pseudo-random number generators. 


The solution of Artin's problem is important for studying the relationship 
between the properties of natural numbers other than +1 and k*2 the 
properties of classes of primes generated by recursive mappings based on 
Fermat's Little Theorem [1,3]. Another relevant applied problem is the 
modeling of self-organizing nonlinear dynamic systems, which are commonly 
called synergetic, taking into account the deep modeling of self-organization 
phenomena in complex systems consisting of sequences of transitions from one 
phase state to another using random number generators with a given 
probability distribution law [3]. 

Numerical sequences of recursive models of cyclic fixed points of dynamical 
systems are determined by the properties of prime numbers, with the help of 
which their behavior is studied under initial conditions 0) =1/ p, where 7 is 
a prime number. In this case, it is necessary to know the law of distribution of 
prime numbers. Riemann proposed the zeta function in 1869: 
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where S is a complex variable, F/ is the set of all prime numbers [1,2]. 
Regarding this Riemann function, a hypothesis was formed, according to which 
all non-trivial zeros of this function are on the line 1/2+ zy, where ;=J/- 1 and 
ye R#., It follows that all primes lie on this line since VY takes values from a set 
that includes all primes 7. Moreover, for any prime number P &(1/2+ ip =O. 
In essence, this was the first attempt to find the law of the distribution of 
prime numbers. 

In 1896, Hadamard and Vallée-Poussin independently proved that the 
equality is true: 
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where (x is the number of primes <-, and the first term in the form of a 
logarithmic smooth function defines the logarithmic distribution of primes in 
asymptotic form however, this information is not sufficient due to the simple 
fact that the logarithmic law does not answer the question: is this number /2 
prime or composite. Enumeration of prime numbers with given properties is an 
even more difficult task [1,2]. In many cases, it becomes necessary to use large 
arrays of primes with given properties. An expressive example of such a set is 
the set of all prime numbers for which a given number 2] is a primitive root, 
provided that this number is not equal to one and is not a perfect square. The 
construction of such sets is a very difficult problem in modern number theory. 

One of the ways to deepen the logarithmic law of distribution of prime 
numbers was the formulation in 1927 by the French mathematician Artin of 
the hypothesis about the primitive roots of prime numbers ye /F, and, 
accordingly, the primitive roots of residue groups (72 modulo a prime 
number / [4,7]. 


Consider the definition of the primitive root of the prime number /. The 
numbers a2+1, a#F is a primitive (primitive) root of the number / if the 
following relations are satisfied: 
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In this case, 22 is a divisor of p- 1=| | 2’. 


7=1 
Given the definition of a primitive root, Artin's conjecture is: 
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where «(x a) is the number of primes / less than or equal to X for which 
a#+l and a+ are according to (1) their primitive roots, da is Artin's 
constant. More precisely, this hypothesis should be presented as follows: 
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But then da, ¥ aa converges in probability to da, and therefore 


has a probabilistic interpretation: da is the probability of choosing from the 
set FP a prime number / such that @ is its primitive (primitive) root. Note 
that the first relation in (4) is always satisfied if @ and V/ are coprime 
numbers according to Fermat's theory [1]. 

It should be noted that Artin offered his estimates for da at a=2. But as 
was proved by Hooley [5], these estimates are not correct. He also proved the 
validity of the relation: 
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with 42) =] [| 1- tot , and the score value is 42) =0.37395581. As will 


be shown later, this estimate is correct only to the accuracy of the first two 
decimal places. 

It should be noted that any number a@a>1 and coprime to V/ is the basis for 
considering the recursive function #¥=a-xmodp, which leads to a recursive 
iterative sequence: 


(6) 


According to Fermat's theorem [1,2], if @ is not a primitive root for /, then 
the process of recursive calculations will continue for such 2/7 that 
Ax, =m =x, -dmodp =1 will be equal ice. 


a” =(modp and z2< p- 1 (7) 


From Fermat's theorem and the properties of the group (7p of residues 
modulo / [1,2], it follows that in this case @ is a generating element of some 
subgroup of the group (24 p4*. Moreover, 7277 is the order of this subgroup, 
which is usually denoted by card p, the number of cosets for this subgroup 
is denoted by “7qd(p. According to the cyclic group theorem (7/ pJ~*, the 
equality always holds: 


p-1 =card p -ind| p (8) 


It follows from the above analysis that Eq. (8) allows us to study Artin's 
hypothesis from a more general point of view, when any natural number a>1 
can be used as a classifier of the set of all prime numbers in terms of ind | D : 
which is the object of further research. As will be established, Artin's primitive 
roots conjecture is a frequent case of its more general formulation. 


II. ANALYSIS OF THE ANALYTIC ASYMPTOTIC METHOD FOR SOLVING THE CLASSICAL 
ARTIN CONJECTURE 


The first attempts to solve Artin's problem are number-theoretic. First of all, 
it should be noted that they are based on the assumption that the generalized 
Riemann hypothesis based on the Dedekind zeta function is true. This does not 
take into account the whole variety of relationships between integers, which 
can be primitive roots, and the set of all prime numbers for which the selected 
integer is a primitive root. The importance of this assumption lies in the fact 
that it is assumed that the Dedekind zeta function takes into account the entire 
set of variants of the dependence of the formation of classes of all prime 
numbers for a given number @ on its properties, which are described by its 
decomposition into prime factors. Although the definition of the Dedekind zeta 
function based on the theory of ideals is based on the Galois theory, this 
dependence cannot be revealed by analytical methods. In addition, only the 
primitive roots of prime numbers are taken into account and the generating 
elements of the subgroups of the cyclic residue group (7 pZ’, as well as their 
properties, are not taken into account. At the same time, information about 
their properties is important, if only because we get new information about the 
structure of the group of residues modulo a prime number VP. 

Probabilistic methods for estimating generalized Artin constants based on 
computer simulation will be considered as another model for solving this 
problem and as a generalization of Artin's original hypothesis. The creation of 
this approach was preceded for a long time by the analysis of the method of its 
solution by the English mathematician H. Holly, who published in the article 


[4], and then re-given for the particular case a=2 in the monograph [5]. First 
of all, let's pay attention to the attitude of well-known specialists in the field of 
number theory to this solution. The Soviet and American mathematician 
Yu. I. Manin in his monograph [2] unequivocally asserts that Artin's conjecture 
has not been proven and is very complicated. The same position is taken by 
such famous mathematicians as Cohen [6], Rosen [7], Moree [8] and others, 
but at the same time C. Pomerance in his monograph [1] notes that this 
problem was solved by Ch. Hooley [4]. Due to this ambiguity of attitude to the 
work of Ch. Hooley and the importance of its solution from both theoretical 
and applied points of view, let us dwell on the works of Ch. Hooley [4] in more 
detail. 

Note that to prove the logarithmic distribution of prime numbers (2), the 
Riemann zeta function (1) was used, but the Riemann hypothesis on the 
distribution of nontrivial zeros of the zeta function was not used in any way. H. 
Hooley in [4] and monograph [5] uses the Dedekind zeta function over Galois 
fields. At the same time, the author assumes that the hypothesis is true: 

Hypothesis: The real part / of every complex zero p=/+7y Dedekind zeta 
function is equal to 1/2 for any Galois field of type G. 

This assumption is equivalent in its complexity to the millennium Riemann 
hypothesis, i.e. Millennium Hypothesis that as provided this assumption is very 
restrictive. The theorem proved by the author has the form: 

Theorem 1. If we assume the validity of the Riemann hypothesis for the zeta 


functions of Galois fields of type OWT, 4/1), where A is without square and 


4‘, | A, then we have: 
a) for any a#0&+1 other than a perfect square, let N(x be the number of 
primes ?<X, for each of which @ is the primitive root of modyp, let also & 


be the without square part of @, let /2 be the largest integer such that @ is 
/2- oh degree of an integer and let: 
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Then, if a £1(mod)), we have: 
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at X> ©, if a =1mod), we have: 
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at X~> ©, 


b) If 2#O,+1 is not a perfect square, then there are an infinite number of 
primes / for which @2 is a primitive root of modyp, 

In the monograph [5] the author repeats the proof only under assumption 
a=2 and the theorem has the following formulation: 

Theorem. If we assume that the extended Riemann hypothesis is valid for 
Dedekind zeta functions over Galois fields and Kummer circular fields of type 
O42, #1), where A is a without square number, then: 


a) Let N(x be the number of primes / not exceeding -X for which 2 isa 
primitive root modulo /. Then: 
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b) There are infinitely many primes / for which 2 is a primitive root 
modulo P. 

In proving this theorem, the author states: 

“Index theory shows that for #2 and any prime divisor @Y of ,~- 1, the 
decidability of the comparison is: 

vu? =Amodyp 

is equivalent to the divisibility of the index of the number 2 (in any base) by 
q7. Thus, we have a criterion that for each 2 is a primitive root modulo / if 
and only if #2, and there is no prime divisor GY of y- 1 for which 2 is the 
residue of the power of GY modulo 2”. 

A comparative analysis of the above two theorems suggests that the author 
has come to the conclusion that the first theorem is not perfect. It is obvious 
that for case a=2 the expression: 
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do not match. According to Hooli's work, / is the maximum positive integer 
for which @ isthe /-th power of an integer modulo P. 

Now let's pay attention to several facts that follow from Hooli's theorems 
and the results of estimating the generalized Artin constants given in the table. 
Hooli's estimate starting from the third decimal place differs from the exact 
estimate obtained from the results of computer simulations. It should be noted 
that the basis of any method for checking whether the number @ is a primitive 
root of a prime number V/ or a generating element of a subgroup of the 
residue group (7p modulo V/ is a recursive function, which is always 
determined by the iterative sequence (6). Ifa natural number @, a#+l anda 
full square, and a prime number V/ are given, then the length of the recursion 
(6) in the general case can only be determined as a result of calculations in 
accordance with this expression. Based on this remark, we can conclude that 
for different values of the number @ for the same prime number /?, the 


lengths of the recursive sequence (6) can differ significantly. The existing 
information technologies for analyzing the behavior of all pairs (a ) are not 
reflected in the generalized Dedekind zeta function. Therefore, an estimate of 


the Artin constant for @ in the case when @=2, obtained by analytical 
methods, can be considered as a value asymptotically close to the exact value, 
but this estimate of the behavior of this constant on the entire set of pairs 
(a D, which is infinite in both @ and P, does not reflect the entire possible 
variety of values. 

From the table below, we can conclude that there is a large variety of Artin 
constants. Hooli's theorems not only do not allow obtaining their estimates, but 
also do not provide information about their values and the dependence of these 
values on the properties of primitive roots determined by decomposition into 
simple factors or their representation in the form @=4+1 or a=4k+3 as 
proposed by Chebyshev [1][2]. Natural numbers other than unity are most 
interesting to consider as classifiers of the set of primes, according to the 
index value determined from expression (8). In this case, the perfect squares of 
natural numbers are the bearer of important information about prime 
numbers. 


III. MODELING THE PROCESSES OF GENERATING DYNAMIC INFORMATION ABOUT THE 
STRUCTURE OF CLASSES OF PRIME NUMBERS ACCORDING TO A GIVEN BASIS 


Now let's return to the logarithmic law of distribution of prime numbers 
[1,2] in order to pay attention to the fact that the given equality does not 
provide comprehensive information about the structure of gaps between prime 
numbers. It is especially important to have information about the distribution 
of smooth primes [1]. This information is especially important when solving the 
discrete logarithm problem and applying algorithms for solving it in modern 
coding theory, modern cryptography. It is known that finding smooth prime 
numbers of large magnitude is very difficult. It follows from this that it is of 
considerable interest to search for the laws of distribution of prime numbers 
not only with respect to their primitive roots, but also with respect to the 
generating elements of subgroups of the residue group modulo a prime 
number (7 p2Z. Artin's hypothesis does not involve such detailed studies. 
Problems of this kind were not considered at all. 

The second circumstance is that simultaneously with this fact, the dynamics 


of change in Quer" is investigated. In [7,8], the entropy of function 


Ax =ARX- Li was estimated and was proved that it has a fractal character. 
These facts are the basis for the formation of proposals on the need to study 
other models of the distribution of prime numbers. 

In addition, it is generally accepted, even at the present time, that it makes 
sense to study it more fundamentally. The first attempt was made by D. Zagier 
[8], but not completed. The results obtained by the author confirm the rather 
complex fractal behavior of this component. Hence it follows that it is 
necessary to significantly improve the study of the depth of classification of 
primes, taking into account all models of the formation of classes of primes for 


any given base a>1. Further more detailed studies of this component confirm 
that although the logarithmic distribution law is fulfilled, nevertheless, 
complete information about the dynamic properties of prime numbers and 
their relationships with their primitive roots still remains poorly studied. In the 
future, we will consider any values of the base @ greater than one. 

According to Artin's hypothesis [4-6], the set of such primes has the 
distribution law (x a) as an expression: 


axa =da x(x (9) 


where (x is the distribution of prime numbers, and da is a constant 
dependent on @. Until now, despite numerous studies, this hypothesis has not 
been resolved. However, it is not known if this is true for any @ values. If the 
hypothesis is correct, then the question remains how to estimate the constant 
da) for each concrete @ and which properties of the number @ influence its 
value. Answers to these questions are still missing. In works [6][7] a detailed 
analysis of all the results of research in the field of solving the Artin hypothesis 
is given. 

It should be noted that the proof of Artin's hypothesis is important both from 
a theoretical point of view in number theory, and from an applied rhenium 
point, because it’s positive solution is important in cryptography, coding 
theory, and the theory of dynamical systems. In [6], a generalized Artin 
hypothesis was formed for any @a>1, i.e. and at the same time @ may not bea 
primitive root. According to Artin’s generalized theory, the following equality 
is true: 


aWxad=dad -nlx (10) 


where a>1, Z is the index of the subgroup of the group (7 pJ of primes in 
the classification of prime numbers generated by the numbers @, da/ isa 
constant. According to the classification built in [6]: 


Kail = pe P\(p-V/cardp =i (11) 


where cara p is the length of the dynamic recursion x, =@2x(modp at 
xX =1, Pis the set of all primes. 
It is not difficult to show that for any a>1 the equality: 


Sv daa =1 (12) 


This means that primes are evenly distributed in classes Aa” for any @. 
By uniformity is meant that within each class of primes Aaz a logarithmic 
law of the distribution of primes is preserved. The constant 4a,/ determines 
the measure of puncturing prime numbers, based on the value @. If 7=1 then 


2 is the primitive root of all primes At a1). For an arbitrary natural number 
* the equality 


aAxad=dain -Alyx (13) 


Moreover, if X> ©, then daz tends to the limit value 427. If we put 
7=1 then da will be Artin's constant for primitive roots. In this case 
az#+l,and a4 fornone ke WN. This is true according to Fermat's theorem 
[1,2]. Wherein, @ is the primitive root of the group of residues (7 p for any 
pe F such that Aa) =| DE P\( D- 0)/ card D =1 . It is important to investigate 
the classes of primes Aa for 7>1 since in this case the positive integer a 
will be the primitive root for the subgroups of the group (7 »Z with the index 
defined by the relations: 


Rail= D\(p- )/card p =ind\ p (14) 


where ind p =/ is the index of the subgroup of (7 pZ. The classes of primes 
Fad have not yet been studied and the distribution of primes in these 
classes is not known. In [1], an assumption was made that Aavz at />1 is 
proportional to AA@ with a factor of 1/z*. Since />1 is considered, in this 
case it is important to know the distribution of prime numbers for the value 
a=k. This is an important generalization of Artin's hypothesis. At the same 
time, the probability of: 


P( pe Rad) = 

pe P\||Rad\/|R=da2 
membership agrees exactly with the provisions of the theory of probability, and 
therefore, estimating 4a on the basis of successive statistical tests and the 
law of large numbers is parity [9-12]. 

The determination of ¢a,/ for any @/ using analytical methods is unlikely 
in the near term. However, the formation and development of experimental 
mathematics [13-15] opens up another way to solve this problem by using 
computer simulation of nonlinear dynamic processes for the formation of 
classes of prime numbers. 

The process of modeling the distribution of primes in_ classes 
Aanv, A1,2),...Aa A,... was reduced to choosing a set of consecutive 
primes from a set of a sufficiently large sample of these classes. The number of 
primes analyzed at each interval of natural numbers was chosen to be 500,000. 
This choice was largely due to the fact that it was previously established that 
reducing this value leads to more significant fluctuations in estimates, 
although convergence to the limit over the entire set of any intervals, even if 
they are not placed consistently, has the same character. 

The process of statistical testing of wef primes for checking their 
belonging to class a7 was reduced to calculating for the selected number 


(15) 


P the recursive procedure ™%=1, %X, =ax(modp until the pairs 
ax =modp were reached at some step 7. Then card p =/ and according 
to Fermat’s lage and ee cyclic group theorem the number 7- 1 is divisible 
by i and then ind(p =(p- )/card p =i, and therefore pe AaJ and if /=1, 
then @ is the primitive root of the cyclic group (7 p27’, and otherwise it is the 
primitive root of some subgroup. At 7>1, we obtain the primitive roots of the 
subgroups of the (7 pZ residue group with the index />1. 

The study of the distribution law of prime numbers / on their belonging to 
Fta,a had the character of consistent statistical tests on the set of natural 
numbers containing the first 500,000 primes. At the first stage, primes / were 
chosen from the set | JAj y oa oe ae With this *= Pooovo. 

For each me 2,....x, we had to solve two problems: check J? for 
simplicity, and if 2= pe /, then Y- 1 was decomposed into simple factors, i.e. 
systematically solved two non-simple problems of checking numbers for 
simplicity and decomposition into simple factors. An effective algorithm for 
solving them was created based on probabilistic methods in the theory of 
elliptic curves. 

As a result of analyzing aq 2,....x, Aal,....A2a/ sets were obtained 
for some /< x and absolutely exact values of their powers were calculated, i.e. 
| A2D|,...| 74a A|, and then estimates of: 


dala =|Aad|/2l2.... dal»x= 
|\AaZ »\|/2( x 


while da1,¥9—> daJ,...dalx— dad with X> © were obtained. 
_ At the next stage, work was also carried out for prime numbers from the 

~ Pooooors++»ZAoo0000 interval and the values of the ¢aJ),...<¢a,/ constants 
were calculated using the same scheme: At the same time / increases. The 
| B>---ABooo000 and | Booooa:- ay ee sequences were combined, and the 
estimates of the generalized Artin constants were again calculated and the 
process of their refinement was studied on the basis of the theory of large 
numbers in probability theory. This procedure continued until 
x= p=1794246' and this is a ten million prime numbers. It was found that 
da,....4a, A in probability converges to some values, the exact values of 
which are irrational and possibly transcendental numbers. In the process of 
estimating the 4a constants, two important theorems were proved: 

Theorem 2. For any aq 2,3.....4... that is not a square, ie. a4 The 
number of non-empty classes of primes tends to infinity at X¥~ ©. 

Theorem 3. For any aq 2,3,..4... that is not a square, ie. azk The 
number of prime numbers in Aa,/ tends to infinity at X¥> ~. 

These theorems are the basis of the convergence of a sequence of statistical 
tests to marginal values. Since for any xe / it is obvious that: 
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at 74, it follows from this that: 
Kk 
>» dad =1 
I=1 


and this is true for all values of X~ %. The review [5] provides an estimate of 
& 2,1), which is identified by 42,1) in our sense, but 42,1) differs from the 
estimate of &2,1) starting from the fifth decimal place and this is a theoretical 
error of the survey works. 

For different 2€| 2,35,6,7,81Q11..|., the behavior of the da constants 
is complex group-theoretic and number-theoretic. The study of their dynamic 
properties is beyond the scope of this work. It should be noted that the results 
of computer simulation of the processes of distribution of primes are 
calculated with an accuracy of the eleventh decimal place for estimates of 
F2,)), d3,0, 65,0, f6,0,.. values. This cannot be asserted for classes by the 
Z7=2 index. To achieve the same accuracy with s=2, it is necessary to 
significantly increase the number of prime numbers. With an increase in the Z 
class index Aaz more than three requirements and the volume of the 
analyzed primes increases in accordance with the unexplored laws. 

Probability-theoretic interpretation of the constant: 


a 


Consider the probability space (2, 7% A based on: 


| | _| a 
C2 =| yp Dare =| Bree Dare =P 


Obviously at X~ © the numbers are z(¥— », a(x a — o, but: 


a(x a =|7a1, | 
a(x =| 709) 


_ |Aal, x| 
cama: 


and at X~ © it is obvious that: 
|Aal, x|/|Ax|— dad 
is where xe P, Po ow, 


Rais = p| ps x&(p-/card p =i 


isat X> © Aaixn— Aad. Thus: 


It follows from Artin's hypothesis that with ca there is precisely the 
probability of a random event -t a1) consisting of a choice of 2 = Q,...,D,-- 
from a prime number / for which @ is an original root of the cyclic group 
(7 pZ . To estimate this probability, the law of large numbers and the method 
of successive statistical tests were used. The essence of the method is that the 
first test group was reduced and calculated for 2D, Diu Dees, by each 
a€|2,3,...16 evaluation of the values of da@zZX at X= Poooo for all 
possible values of *= 12,....4.., that is, Glal4,...Gla A W,... was 
calculated on the next iteration, the same tests were performed for the second 
iteration on the set oer wees Glal, »¥,...6,(41, ¥,... Estimates 


were obtained at the same time G( a1, ¥),....6,(2 & »),..., provided that the 
first and second samples were combined and computed values and were 
determined by |\AazxX- Aal,»|=e for all X. The main focus was on 
c&al,x. As a result of some iterations, it was found that for all @ the 
estimates obtained: 

Ax = pl psx 


Raix= p| p< x&l p-1)/card p =i 


the order of the cyclic group of the subgroup (7p. If /=p-1, then @ is an 
original root, and if /< p- 1 is the original form of the da Artin measure, 
daz isameasure of classes by Aaz in P. At that 627 =|Aa+/|A and at 
the same time: 


x da” =1forall a>1 
i=l 


This applies only to classes with indexes 7=1. For 7>2 it is necessary to 
increase the number of statistical tests. This is naturally due to the fact that 
the classes Aea,z x) for 7>2 from numerical theorems contain less than prime 
numbers. In [1] it is stated that this decrease should be of the order of 1/7 
[15], but this is an erroneous assertion. This is clearly seen from “Table 1”. The 
degree of decline essentially depends on the properties of @ and requires a 
separate study. Case ae 4,916 requires separate investigations, because 
these numbers cannot be primitive roots of that number /, in accordance with 
the Fermat theorem [3] cannot be generating elements of groups (7 pz. 
However, they are generating elements of the subgroups of the group (7 pZ 
with even indices. All classes with odd indices are empty sets. “Table 1” shows 
the constants for &2@1) for all @ except | 449,16. Analysis of the table. The 


table contains over a thousand columns. The analysis of these data is 
numerically theoretical and group-specific and goes beyond the scope. 
The simulation process of the dynamics of the formation of prime numbers 


was constructed on the following assumptions. Suppose that an ordered set of 


; eee : : 
prime numbers P= 2, P,..-.,Q;.--. is given, whose elements are ordered in 


ascending order. All this set was split into a subset of 500,000 primes. The 
number of 500,000 is due to the limitations of MS Excel, as a statistical 
analysis tool, on a number of characteristics of the process of generating prime 
numbers. Only one restriction is important. We always select 500,000 
consecutive primes of the set /. In the current version of Excel, this number 
can be increased to one million. If you use a powerful computer, you can 
choose a larger number instead of a million. 

The implemented version of the study of dynamic processes for the 
formation of primes includes the following indicators: the number of a simple 
number in the / in the ordered set of P, the value of a simple number of /, 
the value of the recursion length of the numbers cara Pp at the same value of 


2 for all prime numbers /P, the index ind\ PD. of the index of the class: 
ind p =| p- )/card p 


the value of the residues modulo any natural module 2>1, for all classes and 
any other analytic properties of primes or factors of the decomposition of the 
number of £- 1 into simple factors. For each simple multiplier 7 in the: 


p- 1=| | 2" 
A 


Decomposition, one parameter of the dynamic process of generating primes 
is presented, with separate indicators that can be analyzed for any other 
indicators, the values for them are deducted by the modulus of the natural 
number z2>1. The only exception is ind\ —. The number of controlled 
indicators analyzed in the Excel environment can be expanded. 

According to the idea of experimental mathematics on the first iteration, 
we proceed from hypothetically known data. But it is also the basis for 
obtaining experimental information on the basis of which the analytical 
methods of the theory of numbers yield an expanded representation of the 
hypothesis in the form /Z,. It is possible that at the same time the hypothesis 
can be corrected or even rejected as not true. From the point of view of 
information technology in mathematics, the hypothesis 7/7, is used to develop 
from the point of view of deepening the experimental mathematics of the 
model of in-depth studies at the level 4. 

The iterations process is continued until an analytically based solution of 
the generated hypothesis is obtained. Since the Artin generalized hypothesis is 
considered in the paper, we present the results of the estimation of the 
constant ¢a@,/ for the case a=4 and /=2. The number a=4 is a perfect 


square, and therefore it cannot be a primitive root. In terms of Artin's 
generalized hypothesis, this is as interesting and important as in the case 
when @ is an original root. 


TABLEI [HE QUANTITY OF PRIME NUMBERS INTO INTERVALS FOR A=2,3,5,8,12 


Interval /a 2 3 5 8 12 
(0.0 - 0.5) 10° 187111 | 187011 | 196980 | 112331 | 187013 
(0.5 - 1.0):10° 186912 | 186948 | 196836 | 112075 | 187057 
(1.0 - 1.5)-10° 186953 | 186960 | 197030 | 112175 | 187040 
(1.5 - 2.0)-10° 186846 | 186856 | 196894 | 112201 | 186958 
(2.0 - 2.5)-10° 187410 | 186896 | 196720 | 112345 | 186792 
(2.5 - 3.0):10° 186711 | 186777 | 196957 | 112042 | 186767 
(3.0 - 3.5)-10° 187096 | 186926 | 197025 | 112335 | 187157 
(3.5 - 4.0):10° 186975 | 187176 | 196942 | 112283 | 186984 
(4.0 - 4.5)-10° 187197 | 187148 | 196543 | 112296 | 187317 
(4.5 - 5.0)-10° 186713 | 186796 | 196689 | 121919 | 186721 
(5.0 - 5.5)-10° 186828 | 187013 | 197050 | 112093 | 187005 
(5.5 - 6.0):10° 187197 | 186771 | 196790 | 112362 | 186936 
(6.0 - 6.5)-10° 186881 | 187116 | 196851 | 112226 | 187056 
(6.5 - 7.0):10° 187065 | 187214 | 196478 | 112093 | 187122 
(7.0 - 7.5):10° 187039 | 186718 | 196957 | 112236 | 187050 
(7.5 - 8.0):10° 187045 | 186756 | 196764 | 112128 | 187161 
(8.0 - 8.5)-10° 187299 | 186805 | 196840 | 112187 | 186594 
(8.5 - 9.0)-10° 186663 | 187050 | 196583 | 111967 | 187144 
(9.0 - 9.5)-10° 186874 | 187156 | 196795 | 112133 | 186976 
(9.5 - 10.0)-10° 187034 | 187072 | 197083 | 111993 | 186947 
10° - 108+0.5:10° | 186932 | 187417 | 196902 | 112048 | 186555 
10° - 10° + 0.5 187123 | 186897 | 197070 | 112386 | 187073 
108 
1012 - 101? + 187127 | 186994 | 196871 | 112303 | 186786 
0.5-10° 


Based on the data presented in [6], we obtained estimates for daz and 
a =2,3,...910and 7=12,..-9.... It is shown that their values are stable for 
class 4,2) i.e. class with iaq(p =2 to within a fourth decimal place. 

The estimates for the 4a constants given in “Table 1” have the unique 
Z=1 property, which is that for aq 235,6,7,101112131415 they 
coincide with the accuracy of the third decimal place. 

The data in “Table 1” allow us to make an important conclusion that there 
are many primitive roots for which the generalized Artin’s constant c& a1) is 
equal to the same value 0.37739... The generalized Artin hypothesis for all 
classes AaW,....Aa7,... will require additional studies based on 
probabilistic computer simulation on the set of prime numbers of data beyond 
the limits of the first hundred million. 

The results of experimental mathematics in “Table 1” of the first iteration 
confirm that Artin's hypothesis is correct. The estimates of the constants are 
obtained with the accuracy of the third decimal place. For 2€ 2,35,6,7,8310 
the: 


» dap =1 
i=l 

and for ae 49 all da2/+D =O and: 
S da22 =1 
f=1 


This is due to the fact that for all @=F this is true because they are 
primitive roots of (7 pZ groups, but primitive roots of their subgroups with 
even indices [3]. 

The results obtained are the basis for constructing an analytical proof of 
Artin’s hypothesis and its general form without using the Dedekind zeta 
function, which will be presented in the next works of the authors. 

The da ratings given in the table for the set of primitive roots 
| 2,3,...16 are obtained for the first time based on the results of computer 
simulation. The literature is known estimation 42,1, which, starting from the 
fourth decimal place, is estimated analytically incorrect, due to the fact that 
the formula: 


1 


&2,1) “IT 1- PGS 


is not true, because it includes all primes and among them those primes for 
which @a=2 is not a primitive root [5]. An important result is the creation of a 
computer model of the process of forming classes Aa@J),....AaJ/,... For 
any values of a>1, the interactions between the classes “Table 2” and 
“Table 3” are investigated (as a continuation). The first estimates were daa 
for 7>2, and it was established that the statement that d4a@/ is proportional 
to 1/7 is absolutely false [1]. Obtaining the results is the basis for further 
deepening research on the Artin's hypothesis using analytical methods. 

The numbers a=3 and a=12 are included in “Table 1” for the reason that 
AZ) and A12D contain 9-8 common primes, while for any other pairs of 
sets A a,,1) and £t a,,1) , the total fraction of common primes is O.4 for /#s 
from the considered set of values. These two facts are obtained on the basis of 
the analysis of the information presented in “Table 3”. This result was also 
based on the methods of modern number theory and probability theory. 


TABLE. TABLE OF GENERALIZED ARTIN CONSTANTS WITH INDEX I VALUES LESS THAN 
OR EQUAL TO 10 


a P(a,1 | P@2 | P@3 | P(@a4 | P(a,5 | P@a6 | P@7 | P@8& | P@9 | P@1 

y) J J J y) y) y) y) y) 0) 

2 0,374 | 0,280 | 0,066 | 0,046 | 0,018 | 0,049 | 0,008 | 0,035 | 0,007 | 0,014 
0 5 4 7 9 8 9 1 4 1 

3 0,373 | 0,299 | 0,066 | 0,056 | 0,019 | 0,033 | 0,008 | 0,014 | 0,007 | 0,015 

9 2 6 1 0 2 9 0 4 1 


4 0 0,560 0 0,093 0 0,099 0 0,070 0 0,028 
9 5 7 1 3 

5 0,393 | 0,265 | 0,070 | 0,066 0 0,047 | 0,009 | 0,016 | 0,007 | 0,028 
7 7 0 4 3 4 6 8 4 

6 0,374 | 0,280 | 0,066 | 0,074 | 0,018 | 0,049 | 0,008 | 0,014 | 0,007 | 0,014 
1 5 5 8 9 8 9 0 4 2 

7 0,374 | 0,282 | 0,066 | 0,068 | 0,018 | 0,050 | 0,008 | 0,017 | 0,007 | 0,014 
1 7 4 4 8 3 9 0 4 3 

8 0,224 | 0,168 | 0,199 | 0,028 | 0,011 | 0,149 | 0,005 | 0,021 | 0,022 | 0,008 
3 3 5 1 4 6 4 1 2 ) 

9 0 0,598 0 0,112 0 0,066 0 0,028 0 0,030 
3 2 6 1 3 

10 0,374 | 0,280 | 0,066 | 0,071 | 0,018 | 0,049 | 0,008 | 0,016 | 0,007 | 0,014 
1 4 5 3 9 9 9 6 4 2 

11 0,374 | 0,281 | 0,066 | 0,069 | 0,018 | 0,050 | 0,008 | 0,017 | 0,007 | 0,014 
1 3 4 5 9 0 9 3 4 2 

12 0,374 | 0,299 | 0,066 | 0,056 | 0,018 | 0,033 | 0,009 | 0,014 | 0,007 | 0,015 
0 1 5 1 9 3 0 0 4 2 

13 0,376 | 0,278 | 0,067 | 0,069 | 0,019 | 0,049 | 0,009 | 0,017 | 0,007 | 0,014 
4 ii 0 Z 1 e) 0 4 4 1 

14 0,373 | 0,280 | 0,066 | 0,070 | 0,018 | 0,049 | 0,008 | 0,017 | 0,007 | 0,014 
9 6 5 7 9 8 9 1 4 1 

15 0,373 | 0,279 | 0,066 | 0,070 | 0,018 | 0,050 | 0,008 | 0,017 | 0,007 | 0,015 
2 6 5 8 9 8 9 7 4 1 

16 0 0,374 0 0,186 0 0,066 0 0,140 0 0,018 
0 9 4 3 9 

17 0,375 | 0,279 | 0,066 | 0,069 | 0,019 | 0,049 | 0,009 | 0,017 | 0,007 | 0,014 
4 4 ie 8 0 7 0 5 5 1 

18 0,374 | 0,280 | 0,066 | 0,046 | 0,018 | 0,049 | 0,008 | 0,035 | 0,007 | 0,014 
0 5 4 7 9 8 9 0 4 2 

19 0,373 | 0,280 | 0,066 | 0,070 | 0,018 | 0,049 | 0,008 | 0,017 | 0,007 | 0,014 
9 8 5 0 9 9 9 5 4 2 

20 0,393 | 0,265 | 0,070 | 0,066 0 0,047 | 0,009 | 0,016 | 0,007 | 0,028 
6 7 0 4 2 4 6 8 4 

21 0,372 | 0,281 | 0,068 | 0,070 | 0,018 | 0,048 | 0,010 | 0,017 | 0,007 | 0,014 
2 9 1 5 8 6 ui 6 6 2 

22 0,374 | 0,280 | 0,066 | 0,070 | 0,018 | 0,049 | 0,008 | 0,017 | 0,007 | 0,014 
0 5 5 4 9 9 9 4 4 1 

23 0,374 | 0,280 | 0,066 | 0,069 | 0,018 | 0,049 | 0,008 | 0,017 | 0,007 | 0,014 
1 8 4 9 9 9 9 5 4 1 

24 0,374 | 0,280 | 0,066 | 0,074 | 0,018 | 0,049 | 0,008 | 0,014 | 0,007 | 0,014 
0 5 5 8 9 8 9 0 4 2 

25 0 0,570 0 0,132 0 0,101 0 0,033 0 0,019 
8 8 ) 3 0 

26 0,374 | 0,280 | 0,066 | 0,070 | 0,018 | 0,049 | 0,009 | 0,017 | 0,007 | 0,014 
1 5 4 2 9 9 0 4 4 Z 

27 0,224 | 0,224 | 0,199 0 0,011 | 0,099 | 0,005 0 0,022 | 0,011 
4 4 4 3 7 4 2 4 

28 0,374 | 0,282 | 0,066 | 0,068 | 0,018 | 0,050 | 0,009 | 0,017 | 0,007 | 0,014 
0 8 5 4 8 3 0 1 4 3 

29 0,374 | 0,280 | 0,066 | 0,070 | 0,018 | 0,049 | 0,008 | 0,017 | 0,007 | 0,014 
5 1 6 0 9 8 9 6 4 1 

30 0,374 | 0,280 | 0,066 | 0,069 | 0,018 | 0,049 | 0,008 | 0,017 | 0,007 | 0,014 
0 5 5 9 9 9 9 8 4 y 

31 0,374 | 0,280 | 0,066 | 0,070 | 0,018 | 0,049 | 0,008 | 0,017 | 0,007 | 0,014 


1 6 5 1 8 9 9 5 4 2 
32 0,295 | 0,221 | 0,052 | 0,036 | 0,094 | 0,039 | 0,007 | 0,027 | 0,005 | 0,070 


3 4 4 9 ) 4 0 7 8 9 
53 0,374 | 0,280 | 0,066 | 0,070 | 0,019 | 0,049 | 0,009 | 0,017 | 0,007 | 0,014 
0 4 5 1 0 8 0 5 4 2 


In conclusion, by returning attention to “Table 1” of another theory of vision. 
The essence of a fundamentally new fact is that wherever 500,000 primes 
be F are selected for any @=>, the number of primes in classes ranges from 
no more than 500, which is no more than a thousandth of them. This means 
that on any set of consecutive primes we obtain an estimate of the Artin’s 
constants up to the fifth decimal place. Statistical summation of values over 
the entire set of the first ten million primes made it possible to obtain 
estimates of the constants c& a1) to the eighth decimal place. 

It follows that the methods of computer simulation of the processes of 
forming classes of primes Atal), Aa@2),...Aaad,... and estimation of 
constants &4al),da@,2),...daJ/,.. are the basis for the development of 
information technologies in modern both pure and applied mathematicians. 

“Table 2” shows the values of the estimates of the generalized Artin’s 
constants for the marked set of values of a, which were studied as classifiers of 
the set of all primes. 

“Table 2” shows a fragment of the entire huge number of obtained estimates 
of the Artin’s constants. The first column contains estimates of Artin’s 
constants in its original form. They belong to primitive roots. The numbers 
4,9,1 6 are not primitive roots, since according to Fermat’s theorem they, like 
the squares of other numbers, cannot be primitive roots. However, they can be 
classifiers of primes as roots of subgroups of residues modulo primes. An 
interesting fact is that they can be used to build pseudo-random number 
generators. In addition, the diskette logarithm problem can be considered on 
their basis. 


IV. DYNAMIC PROPERTIES OF THE FORMATION OF CLASSES OF PRIMES IN THE 
GENERALIZED ARTIN CONJECTURE 


In accordance with the developed mathematical model for the formation of 
classes of primes by base a@>1 and the calculated values of the generalized 
constants ¢a,/ for 7=1, as a result of computer simulation, it was found that 
the generalized hypothesis is correct. “Table 1” shows the values of Artin's 
constants, relationships between classes, the dynamics of class formation and 
its properties on the set of all primes P. 

In reality, the modeling of the processes of formation of classes Aa was 
carried out for the set aq 2,3,...172953. The numbers ae 4916 as 
squares of numbers, according to Fermat's theorem [2], cannot be primitive 
roots of p<, and, accordingly, the residue groups of (7p modulo LP. 
Particular attention was paid to the numbers | 5,131 72953 due to the fact 
that they belong to the class of numbers of the Chebyshev type [1,2], that is, 
they have representations »p=44+1, while pe /, and the number J7 is a 
natural number. According to Chebyshev's assumption, the behavior of these 


numbers in the classes of residues modulo a prime number must differ from 
other prime numbers. 

To solve the problem of modeling classes of primes on a given basis and 
estimating the generalized Artin constants daa, an Excel-based software 
package was created that allows you to extend the modeling process to any 
natural numbers a>1, and any set of consecutive primes whose power is a 
multiple of 500,000. This is the number of primes numbers was chosen 
because it is_ statistically representative and provides an _ accurate 
representation of the dynamic processes that form the Aa classes. “Table 
1” shows the results of the simulation process for 2a€ 2,35,812 values. This 
set includes a=2 for the reason that it can be seen that the estimate [5,6] 
differs from the exact value. The difference starts at the third decimal place. 
This fact is important due to the fact that expression (4), although from an 
asymptotic point of view, is close to the exact value of 4&2), nevertheless, it 
does not take into account all the features of the formation of classes Ata) 
for a=2. The number @=5 is interesting because a2=5=4-1+1 is the smallest 
Chebyshev number that is most sensitive to the established fact that all classes 
F051 0k+ 5) for k=O are empty. This is true for all Chebyshev numbers. The 
proof of this fact is number-theoretic, and therefore, is not given. 

The numbers 82'732 are interesting because the dynamical properties of 
the classes A8,2 are radically different from those of the other classes 
studied. In particular, it is established that if a=8 is a primitive root of pe * 
then a=2 is also a primitive root of the same prime number. Conversely, if 
a=2 is a primitive root of ye /, then a=8 will be either a similarly primitive 
root of 2 or pec A83). This is brand new information about Artin's 
generalized constants. The developed approach made it possible to obtain 
fundamentally new results in modern number theory, and, as a result, in 
modern cryptography. 

“Table 2” shows estimates of the generalized Artin constants. The first 
column contains a list of 2& 2,3,...3253 values, as an indication that this is 
reaching the limit of prime numbers of type =44+1. The second column and 
subsequent columns correspond to Artin's constants in his interpretation for 
the z2da, P< 1,2,...10 index values. The second column gives estimates of 


Artin's constants for values of @ as a primitive prime root of ten million 
consecutive primes. Starting from the third column, estimates of the 
generalized Artin constants are found for the reduced index values equal to the 
number 2 and up to 1C. A preliminary analysis for the numbers @a=5 and 
a=2C shows that the dynamics of all generalized Artin constants is almost 
identical. The existence of windows with zeros for these values is a simple 
rationale for increasing the values of Artin's constants, which was not justified 
at the time of the analysis of this problem by both Artin and Bilhartz. 
Moreover, at @=8 and a=27 the Artin’s constants are the same, but the 
behavior of 72da, p >1 these values is very different. 

Let's take a look at @=32. The generalized Artin constants take on 
completely different values. This implies that the theory of the Artin hypothesis 


in a generalized form requires a significant deepening and development for all 
values of a> 2. 

In closing, looking back at “Table 1” from a different perspective. The 
essence of a fundamentally new fact is that no matter where 500,000 prime 
numbers ye Ff are sequentially allocated, for any @ more, the number of 
prime numbers in classes varies within no more than 500, which is no more 
than a thousandth of them. This means that on any set of consecutive primes 
we obtain an estimate of Artin's constants up to the fifth decimal place. 
Summation by statistical methods of the values over the entire set of the first 
ten million primes made it possible to obtain estimates of the constants & a) 
with an accuracy of eight decimal places. 

Hence it follows that the methods of computer simulation of the processes of 
forming classes of prime numbers Aa), Aa2),...Aad,.. and estimating 
constants 4al),da@2),...da/,.. are the basis for the development of 
information technologies in modern both pure and applied mathematics. 

An interesting result is that the constants: 


E2,) =h3,)D =d6) =47.D) = 
d1 QD =...=d15)D =d1 7D... 


are equal to one thousandth, although 48,1 and 45,0 are radically different 
from them. On the basis of modern number theory and the theory of random 
processes, the validity of such results is proved. The proofs of these statements 
can be built only on the basis of data obtained as a result of computer 
simulation. 

Now consider for all prime numbers of the given sets the dynamics of the 
change in the number of different prime factors into which the numbers p- 1 
are decomposed. Work in the Excel environment when creating a modeling 
environment in which all information is available regarding the set of all 
characteristics of the processes of formation of generalized Artin constants. 
“Table 3” shows the dynamics of changes in the value of «(p- 1). In the first 
three columns for  p- 1), the numbers of primes decrease, but starting from 
cA p- ) >4, these numbers increase. You should pay attention to the simple 
fact that the more o{p- ) the greater the rate of increase. But this means 
that the number of primitive roots in sets of consecutive primes will decrease. 
If this assumption is correct, then Artin's hypothesis will have to be 
substantially corrected. 


TABLE Ill. DISTRIBUTION OF VALUES a p- D 


erage w(p- | 5 3 A 5 6 7 | 8 9 
(0.0 - 0.5) 10° |5038 1676 |18870 | 8485 11143]299 |- : 
1 91 1 1 7 


(0.5- 1.0) -10° | 4485 | 1535 | 18825 |9485 |1765|772 |2 - 
2 74 7 6 7 


(1.0-1.5) -10° | 4314 | 1496 | 18688 | 9890 | 2032 | 1125 |6 - 
3 21 0 1 A 
(1.5 - 2.0) -10° | 4204 | 1470 | 18600] 1014 | 2207 | 1374 | 8 - 
9 o1 6 36 6 
(2.0 - 2.5) -10® | 4133 | 1453 | 18546 | 1029 | 2337] 1517 | 14 - 
0 33 0 72 4 
(2.5 - 3.0) 10° | 4078 | 1436 | 18524 | 1044 | 2412 | 1736 | 16 - 
0 01 ©) 99 3 
(3.0 - 3.5) 10° | 4032 | 1426 | 18437 | 1058 | 2492 | 1926 | 20 - 
6 86 2 23 7 
(3.5 - 4.0) -10° | 4016 | 1413 | 18417] 1065 | 2578 | 2028 | 36 - 
.) 04 Z 08 ye. 
(4.0 - 4.5) -10° | 3969 | 1408 | 18381 | 1070 | 2632 | 2187 | 28 - 
4 83 1 71 6 
(4.5 - 5.0) 10° | 3923 | 1398 | 18419] 1073 | 2687 | 2276 | 35 - 
il 96 2 o1 9 
(5.0 -5.5) -10® | 3904 | 1396 | 18348 | 1081 | 2723 | 2349 | 36 - 
3 88 4 62 8 
(5.5 - 6.0) -10° | 3896 | 1389 | 18324 | 1084 | 2792 | 2487 | 37 - 
0 09 7 35 .) 
(6.0 - 6.5) -10° | 3862 | 1384 | 18276 | 1092 | 2834 | 2563 | 51 - 
6 12 8 37 3 
(6.5 - 7.0) 10° | 3851 | 1375 | 18279] 1098 | 2868 | 2560 | 52 - 
3 93 1 04 7 
(7.0 - 7.5) -10® | 3854 | 1367 | 18283] 1099 | 2913 | 2705 | 42 - 
4 88 6 995 0 
(7.5 - 8.0) -10® | 3801 | 1368 | 182331106 | 2937 | 2759 | 55 - 
1 99 3 12 1 
(8.0 - 8.5) -10® | 3799 | 1360 | 18220]1108 | 2947 | 2853 | 62 - 
1 71 9 03 .) 
(8.5 - 9.0) -10° | 3781 | 1362 | 18218)|1109 | 2982 | 2887 | 57 - 
4 40 6 93 3 
(9.0- 9.5) -10® | 3798 | 1354 | 18214 ]1111 | 3019 | 2952 | 66 - 
8 fas) 2 87 2 
(9.5 - 10.0) -10°® | 3733 | 1356 | 18194 )|1117 | 3025 | 3003 | 58 - 
0 88 9 20 2 
10®°- 108+0.5 |) 3305 |1225 |17597)]1211 | 4071 | 6235 | 328 | 2 
10° 2 84 0 19 0 
10°-10°+ 0.5 | 2945 | 1121 |17028)1275 | 4990 | 9791 | 810 | 22 
10° 3 82 4 o1 d 
10-10% 4+ 0.5 | 2233 | 9132 | 15322 | 1363 | 7017 | 2157 | 364 | 348 
10° 9 4 2 64 6 .) 3 


wl p- D is the number of different prime divisors of the number Z- 1. 
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Based on the analysis of the processes of formation of classes of prime 
numbers, fundamentally new information technologies for solving complex 
mathematical problems using the methods of modern’ experimental 
mathematics were created for any bases. The correctness of the developed 
approach and computational efficiency are proved. A generalized theory of the 
Artin hypothesis has been developed, in which its classical version is a very 
special case. Estimates of Artin's constants for bases greater than the number 
two are obtained, and the statistical consistency of the obtained estimates is 
proved. A detailed analysis of the classes of prime numbers has been carried 
out and the foundations of effective methods for the structural analysis of 
classes have been created. It is proved that a new method for modeling the 
dynamics of the formation of classes of prime numbers and describing their 
properties creates the basis for constructing more advanced models of pseudo 
prime number generators, the development of new methods for protecting 
information in modern cryptography, and opens up new opportunities for 
building models of nonlinear dynamic systems. 
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